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• Crystals of He contain vacancies that move around by a quantum mechanical hopping process. 

^\ I The density and pressure of these vacancies can be experimentally studied. The accuracy of the 

experiments is high enough to detect the effect of the Bose statistics of the vacancies. In this paper 
we examine the effect of the hard-core repulsion between the vacancies, which should also have a 
measurable effect on their behaviour. We set up a virial expansion for a lattice gas of hard-core 
■ particles, and calculate the second virial coefficient. It turns out that the vacancies behave as ideal 

Bose particles at low temperatures, but that the hard-core interaction makes them behave more and 
more like fermions as the temperature increases. 

PACS numbers: 67.80.Mg, 05.30.Jp, 61.70.Bv 
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I. INTRODUCTION 



m . 

^\ , The thermodynamic behaviour of vacancies in solid 
*He is an interesting experimental and theoretical prob- 
lem [0. Vacancies in Helium are more mobile than in 
^ ' any other solid. At the temperatures where sohd Helium 
exists the motion of vacancies requires a quantum me- 
chanical description. They hop from site to site with a 
Q certain rate Vy, leading to a band of states e(fc), much 
O hke the electron motion in the tight binding approxima- 
^ tion. For "'He the vacancies are obviously bosons, since 
• i-H , the creation operator for a vacancy is the annihilation 
' operator for a '^He particle, which is a boson The 
^ , fact that two vacancies cannot occupy the same lattice 
■ - - ' site has to be incorporated as a hard-core potential for 
the hopping bosons. So the simplest model for vacan- 
cies in ^He is that of a gas of hard-core bosons on a 
lattice. In reality the strain fields around the vacancies 
produce a more complicated interaction between them 
than the simple on-site exclusion. However, we consider 
the hard-core boson approximation as a sufficiently real- 
istic description of the vacancy motion to leave out these 
further refinements, in order not to complicate the model 
too much. The hard-core bosons on a lattice are known to 
be equivalent to a spin-i quantum system with an inter- 
action of the XY type Q . It is however not this analogy 
that is exploited in this paper. The reason is that the 
vacancy system is in practical circumstances always ex- 
tremely dilute, which corresponds in the spin analogy to 
a system in an extremely large external magnetic field. 
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In this limit the spin analogy is of little use and in our 
opinion less transparent than the particle language. 

Experiments show that the percentage of vacancies in 
a crystal is at most of the order of 1% Q], at least at 
temperatures of the order of IK, where the experiments 
take place. As a first approximation the vacancies thus 
behave as an ideal gas. However, present day experiments 
are sufficiently accurate that effects of Bose statistics can 
be detected. It is one of our points that then also effects 
of the hard-core interaction become detectable. 

In this paper we present a systematic analysis of the 
vacancies using exp(— A/fc^r) as a small parameter, A 
being the excitation energy, or band gap, required to cre- 
ate a vacancy. This is equivalent to a virial expansion |^] 
for the quantum lattice gas, and we work out the proper- 
ties in detail up to the second virial coefficient. A general 
formula for the second virial coefhcient of the hard-core 
Bose lattice gas is derived and evaluated for hypercubic 
lattices in one, two, and three dimensions. The case of 
the hep lattice, which applies to real solid "^He, will be 
treated elsewhere. 



II. THE HARD-CORE BOSE LATTICE GAS 

The vacancies are represented by Bose creation and 
annihilation operators b] , bi obeying the usual Bose com- 
mutation relations. The Hamiltonian for the vacancies is 
given by 



H = -tY^iblb, + blbi) + {A + ct) bib, + - ^ blblbA. 



Here the transfer integral t is equal to ft.z/„, and the hops 
take place between all pairs of nearest neighbours 
on the lattice. The coordination number of the lattice 
is c, and A is the energy required to create a vacancy. 
It functions as minus the chemical potential for the va- 
cancies. The last term represents the vacancy- vacancy 
repulsion. It could be omitted in favour of a change to 
on-site Fermi commutation relations for the b[^\ but we 
prefer to work with standard commutation relations and 
a potential U which penalizes the simultaneous occur- 
rence of two or more vacancies at the same site. We will 
let t/ — !■ cx), or [/ is much larger than any other energy 
in the problem. 

Without the potential term we have an ideal Bose lat- 
tice gas, the Hamiltonian of which can be diagonalized 
by the canonical transformation 



(2.1) 




(2.2) 
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where N is the number of sites in the system. We use 
periodic boundary conditions, and the wave vectors k are 
restricted to the first BriUouin zone of the lattice. 



Inserting (??) into the first two terms of (2.1) we obtain 
the unperturbed Hamiltonian 

Ho=Y.[A + e{k)]b^{k)b{k), (2.3) 



where e{k) is given by 

e{k) =t^{l- cos k-fs), (2.4) 
s 

and fg is the set of the c nearest neighbour positions 
with respect to a centrally chosen site. A is the gap of 
the energy band, since we have £(0) =0 as lowest energy 
in the center of the BriUouin zone. One sees that t is 
a measure for the bandwidth, which for a d-dimensional 
hypercubic lattice is -u; = Mt. 

The interaction is written in terms of 6*^^^ (fc) as 

^ - ^ E hM.+>^.+GbHkl)bHk2)b{h)b{k,)^ 



2N 

ki,G 



(2.5) 



where G is a vector of the reciprocal lattice. The ma- 
trix elements of the interaction have no other structure 
than the conservation of the total incoming and outgoing 
momentum (up to a reciprocal lattice vector), a feature 
which is of great advantage in solving the two-particle 
problem. 



III. THE VIRIAL EXPANSION 

The grand partition function of the vacancy system is 
given by 

S = Tre-^^, (3.1) 



with Ti, given by (2.1), (3 = l/fc^T, and Tr stands for the 
trace over all symmetrized states. As 7i conserves the 
number of vacancies, S can be expanded as 



S = ^ Z„e-"^^, (3.2) 



ri=0 



where Z„ is the canonical partition sum for n vacancies 
excluding the contribution from the gap A. Of course, 
Zq = 1, and 

Zi = Tn e-^^^"-'^) = ^''''^'^ ' (3.3) 



because for one vacancy no hard-core effects enter. Tr„ 
is the trace over n- vacancy states. 
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For In,:^ we may deduce from (3.2) 

oo 

PpNvo = InS = N^bee-'^'^^, 



(3.4) 



where vq is the volume of the unit cell, p is the pressure 
of the vacancies, and the be are the fugacity expansion 
coefficients. The first of these, hi, reads 



(3.5) 



BZ 



where we have replaced the sum over k by an integral over 
the Brillouin zone. In contrast to the continuum ideal 
gas, bi is not simply related to a thermal wavelength A 



(3.6) 



A formula of this type only results when the temperature 
is so low that e{k) may be replaced by its low-momentum 
behaviour 



e{k)^'-J2(^-fsr^^., 
2 ^ 2m* 



(3.7) 



defining an effective mass m* for the vacancies. For d- 
dimensional hypercubic lattices we find 



m 



(3.8) 



where a is the lattice constant. Using (3/7) in ( |3.5[ ) and 
extending the /c-integral beyond the Brillouin zone t o in - 
finity (as is allowed for large /3 or small T) one finds ( |3.6| ) 
with 



= h^/2Trm*kBT. 



(3.9) 



The second term 62 in (|3.4|) is our main concern in this 
paper. We find from (3.2) and (3.4) 



b2 = -(z2-^ 

N \ 2 



l|Tr2e-''(«-^)-i(Tri, 



(3.10) 



We rewrite this expression by adding and subtracting the 
contribution of the non-interacting two-vacancy system. 
So we define 



g-/3CH-2A) _ g-/3CHo-2A) 



(3.11) 



as the contribution of the hard-core interaction, and 



N 



Tr2 e-^(^«-2A) „ 1 l^Tn , 



-/3CHo-A)^ 



(3.12) 
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as the effect of tlie Bose statistics of tlie vacancies. Tlie 
combination 



tint 



(3.13) 



yields the total effect. 

The statistical effects are trivial to calculate, as the 
unperturbed grand partition function is given by 



lnS" = -;^ln(l- 



such that 



(3.14) 



(3.15) 



with &i(/3) given by ( |3.5[) . Thus it suffices to focus our 



attention on the calculation of b™'^ . As a gene ral o bser- 



vation we note that the two contributions in ( 3.13 ) will 
have opposite signs. To see this more clearly, we go over 
to a series in the density n of the vacancies, 



1 1 °° 



-il3A 



=1 



(3.16) 



Eliminating e ^'^ from (3.4) and ( 3.16| ) we obtain a virial 
expansion for the pressure 



(3.17) 



So at fixed density n the statistical effect s lower the pres- 
sure as &2 is positive according to ( 3.15 ). The hard-core 
repulsion can only increase the pressure, so 6™* must be 
negative. 



IV. THE SECOND VIRIAL COEFFICIENT 



The second virial coefficient ( 3.11 ) is evaluated as 
where the level densities P2{E) are given by 



P2{E)=i:v2 5{E-n) 
pl{E)=Tv2 5{E-no). 

The level densities are obtained from the formula 



(4.2) 



Q+{E) 



1 



E + ie-H 



V- 



1 



E-n 



i7rS{E-H), (4.3) 



which transfers the problem to the determination of the 
Green's function, for which the general equation holds 
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giz) = go{z) + go{z)vg{z), 



(4.4) 



with z a complex number. The density p{E) is obtained 
from g as 

p{E) = --ImTrg+iE). (4.5) 

TT 

The states of the unperturbed 2-vacancy Hamihonian 
are denoted by two wavenumbers ki and fc2 with the 
property 

Holh k2) = £'o(fci,fc2)|fci ^2) 

= (2A + e(fci)+£(fc2)) |fcifc2). (4.6) 



The matrix elements of V are obtained from (2.5) 

(feifc2i^ifc;fc^) = ^E%.+fe.,fe;+^+G' (4-7) 

G 

which shows that a representation in center of mass and 
relative coordinates will be advantageous. Thus we in- 
troduce 

K = ki+k2, k = {ki-k2)/2. (4.8) 

The sum over G is eliminated by choosing the Brillouin 
zone in such a way that no two points in it have values 
of K differing by a reciprocal lattice vector, so that only 
the term with G — contributes. The matrix element of 
V is then diagonal in K. From now on we will assume 
that the Brillouin zone has been chosen in such a way, 
and drop the reference to G. Since the total momentum 
K is conserved by Hq and V, g becomes diagonal in it. 
So, writing (4.4) in the K,k representation, we have 

{Kk\g(z)\Kk') = ^ ^ ^ I Srp 

^ ^ ^' ' z-Eo{K,k) \ 

Of/ 

+ ^J2(^k"\g{z)\Kk') 



N 

k" 



(4.9) 



The simplifying feature of (4.9) is that the general matrix 
element of g couples only to the total sum over the first 
entry of the matrix elements. For the latter we obtain an 
expression by summing (4.9) over k 

y^{Kk\g{z)\Kk') 

Y z-E^{K,k') 

+2f/7^(z, K)^{K k"\g{z)\Kk'), 

k" 

(4.10) 

with TZ{z, K) given by 
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n{z,K) 



-E 



- z - 



1 



(4.11) 



Now pl0| ) is an algebraic equation for the quantity in 
t he l eft hand side. Using the solution of this equation in 
(4.9) one finds 



{Kk\g{z)\Kk') 



1 



z-Eo{K,k) 
2U/N 1 



k.k' 



1 - 2UTZ{z,K) z- Eo{K, k') 



(4.12) 



Note that the off-diagonal elements are of order iV~ ^ 
while the diagonal elements are of order 1. With ( 4.12 ) 
we can calculate the level density from 



Trg(2) =^(X^|g(z)|i?fc), (4.13) 

K,k 

leading to the compact expression 

TTgiz)^TTgoiz) + -^J2^ii[l-2UTZ{z,K) 

K 



(4.14) 



which in turn yields for the difference p2 — p\ 

P2{E)-pI{E) = 
d 



1 

— Im 

TT 



K 



z=E+ie 

(4.15) 



This expression holds generally for any on-site repul- 
sion U. We may let [/ — > cx3, by which it will disap- 
pear from the formula as the U term under the logarithm 
starts to dominate the argument for any z and K. Omit- 
ting the 1 in the argument of the logarithm, the term 
ln(— 2f7 ) dro ps out after differentiation with respect to z. 
So for (4.15) we have in the limit U oo the equivalent 
expression 



1 



P2{E)~p',{E) = —lm 



d 



Y.\nniz,K) 



K 



(4.16) 



Hereby the problem is essentially reduced to the evalua- 



tion of 7^(z, K) given by ([l.tl]) . 

A few comments are in order about this expression. 
For bosons the state |fci ^2) is the same as 1^2 fci). So 
the relative momenta k and —k should be identified with 
each other. Both ki and k2 run through a Brillouin zone 
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appropriate for the structure of the lattice, and this in 
principle defines the ranges of K and k. But as was 
mentioned before, we have chosen the Brillouin zone such 
that there are no points whose values of K differ by a 
reciprocal lattice vector. 



From the definition (4.11) of TZ{z,K) and the expres- 



sion (4.6) for Eq{K, k) one sees that 2 A occurs as a shift 
in the energy variable and the p2 are functions of E—2/S.. 
Taking E—2/S. as an integration variable in (4.1) one finds 



1 

77 



dE' [p^iE') - pI{E')] 



-I3E' 



(4.17) 



with 



P2{E') - pI{E') = 

TT 



and 7?.'(z, K) given by 



K 



n'{z,K) 



e(f + fc)-£(f -fc) 



(4.18) 



(4.19) 



In this formula A is eliminated, as it should be. U sing t he 
fact that Imlnz ~ argz, we can further rewrite (4.18) 



P,{E')^pI{E')^--^F{z)\ 

7T OZ 



z — E'-\-ie 1 



(4.20) 



with 



K 



(4.21) 



From (H.18|) and ( |4.19| ) one sees that only E' values 
occur in ( 4.17 ) which lead to complex values of TV{z, K). 
These occur when z = E'+ie is a pole in the /c-integration 
in ( [4.19| ). Thus the combined bandwidth of e{K/2 + k) + 
e{K /2 — k) determines the range of E' values. This is 
twice the ba ndwid th of £{k) just as in the ideal Bose 
contribution ( ^.15| ). 



V. THE ONE-DIMENSIONAL CASE 



We interrupt the general disc ussio n for t he treatment 
of the one-dimensional case of (4.17)-(4.21 ), as this case 
is interesting, completely analyzable, and elucidating for 
the structure of the functions occurring in (4.17)-(4.21). 

The d = 1 band structure is given by (for a lattice 
constant a = 1) 



e{k) = 2t(l - cos A:), 
which gives for 6i(/?) the expression 



(5.1) 
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6i(/3) = r dfce-2/3t(i-cosfe) ^ 5-2/9*70 (2/3t), (5.2) 

271" J^n 

where /o(a;) is a modified Bessel function of tlie first kind. 

For the evaluation of TZ{z,K) we rearrange the Bril- 
louin zone in such a way that its boundaries in K, k space 
are convenient, and also such that no reciprocal lattice 
vector G enters into the problem. In figure we have 
divided the original Brillouin zone — tt < fci < tt , — tt < 
k2 < TT into four domains, I, II, III, and IV. Domains I 
and IV, as well as II and III, refer to the same states, as 
they are obtained from each other through interchanging 
fci ^ fc2. So it suffices to take one of each, say I and III. 
Now III is equivalent to III', which follows from III by 
shifting fci over 27r. The combined domain I and III' is 
given in K, k space by0<iir<27r, 0<fc < tt. 

Using this parameter space one can write (4.19) explic- 
itly as (dropping the primes) 



1 r 
n(z,K) = — / dk- 



1 



4i(l — cos Y cosk) 
For convenience we put z = 4t(C + 1 + ie), and have 

I r „ 1 



(5.3) 



SlTt 



dk 



( + ie + cos Y cos k 



(5.4) 



According to equations ( 4.20| ) and ( 4.21 ) we are inter- 
ested in real values of (. So, letting e i 0, the integral 
has a real part given by a principal value integral, and an 
imaginary part given by an integral over a delta function 



C -I- cos Y cos k 

i K 
— I dk S{( + cos — cosk). 
8t Jq 2 



(5.5) 



Two cases must be distinguished: for |C| > | cos(iir/2)|, 
the imaginary part of TZ is zero, and for |C| < | cos(i4r/2)| 
its real part is zero. The result is 



1 



sgn(C) 



8t 



„2 K 



(|C|>|cos(if/2)|) 



8t 



cos^ — C"^ 



(|C|<|cos(if/2)|). 



(5.6) 



The important quantity in equations ( 4.20 ) and ( 4.21 ) is 
arg7^. From (|!|) we see that if C < -| cos(is:/2)|, 7^ is 
real and negative, while for C > | cos(if/2)| it is real and 
positive. For -| cos(if/2)| < C < |cos(if/2)|, 7^ is pure 
imaginary with a negative imaginary part. So we find 
that 



<i>(C,i^) = arg7^(C,if) = 



(C<-|cos(if/2)|) 
-n/2 (Id < |cos(if/2)|) 
(C> |cos(if/2)|). 



(5.7) 
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These values of $ are plotted in figure || in the C,,K/2 
plane. Integrating in the JC-direction, we find for -F'(C) 



K "^0 



-TT 

- arccos ^ 







(C < -1) 

(-i<C<i) 

(C>i). 



(5.8) 



Using C instead of £' in (4.17) as integration variable we 
have 



^2'^^ f^dC [P2(C)-P§(C)] e- 



-4/3t(C+l) 



(5.9) 



Due to the differentiation with respect to z in ( 4.20 ) only 
values of ( between —1 and 1 contribute to the integral. 
From (5.8) we see that for |C| < 1 



P2(C)-P^(C) 



N d 



F{0 



N 



Substituting this into (5.9) we find 
= -e-4/5*/o(4^t). 



(5.10) 



(5.11) 



We now compare this with the result for an ideal lattice 
gas, given by ( [3.15D 



lilt 
'2 ■ 



(5.12) 



So the hard-core interaction giving rise to 62"* changes the 
value of 62 from into its opposite, or in other words, the 
Bose value is turned into the Fermi value. This is exactly 
what has to be expected from the well known fact that a 
hard-core Bose gas in one dimension is equivalent to an 
ideal Fermi gas (for all virial coefficients). 

The picture shown in figure ^ for the phase $ has some 
general validity, in the sense that for sufficiently negative 
C one has <I> = — tt while for ^ sufficiently positive $ = 0. 
In the zone in between one has for d > 1 in general a 
continuous transition from — tt to 0, with eventually also 
zones with $ = — 7r/2. 



VI. THE TWO- AND THREE-DIMENSIONAL 
LATTICES 

The calculation of 62 for hypercubic lattices in higher 
dimensions runs along the same lines as what was done 
in section for one dimension. The band structure is 
now given by (again a = 1) 

,7^, _ J 2t(2 — cosfc^: — cosfcj,) {d — 2) 

^ ' 1 2t(3 — cos fcj: — cos fcj, — cos fcz) {d — i). 

(6.1) 
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In d dimensions, 6i(/3) is simply given by 



(27r)" 



^;^g-2/3t(l-cosfc) 



= e-^'^'^' [Io{2(3t)f . (6.2) 

As before, we rearrange the Brillouin zone such that its 
boundaries in K, k space are convenient. The different 
components of K, k are independent, and we have < 
Ki < 211 ^ < kx < TT , — TT < ky,kz < TT. The kx- 
interval is halved to avoid counting the same (symmetric) 
state twice. The two-vacancy energy bands are given by 

Eo{K, k) = 4t(2 — cos cos k^ — cos cos ky) 

(d=2), 

Eo{K, k) = 4t(3 — cos cos kx — cos cos ky 
— cos cos/cz) {d — 3). 



(6.3) 



We scale the parameter z as 



8t{C + l + ie) (d = 2) 
12i(C+l + ie) (d = 3). 



(6.4) 



so that in both cases the energy band runs from = — 1 
to C = 1. 



In two dimensions wc find with (5.5) 

1 



(2^)2 



dk^ 



1 



1 



32t7r 



8i 
f{A,B) 



Kx 



C + ie + - ( cos cos kx + cos cos ky 



(6.5) 



where 



1 K 

A = C + 2 



1 if. 
- cos — — 

2 2 



(6.6) 



and 



f{A,B) 



sgn(A) (A2 > S2) 
-^ (^2 < S2). 



(6.7) 



The integral in (3.£) can be expressed in terms of com- 
plete elliptic integrals of the first kind (see appendix). 
This gives an analytic expression for 7?.2(Cj K). 

This expression can also be used to find the result in 
three dimensions, 



7^3(^,i?) = 



1 



(27r)3 



dk. 



dk. 



Kx K Kz ' 

4t(3 — cos cos kx — cos cos ky — cos cos k^) 



(6.8) 
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dk. 



ZTT 



C + + - (cos cos kx + cos cos fcy ) 



(6.9) 



where 



8t{C + 1) = 12t(C + 1) - 4t(l - cos -y COS fc^). (6.10) 

Using the analytic expression found for TZ2, TZ^ can be 
calculated by numerical integration of (6.9). 



The next step is to obtain F{(^) as given in ( 4.21 ) by 
integrating over K. This is done numerically using Monte 
Carlo integration. The result is a function F{C) that is 
equal to — vr for C < — 1, where TZd is real and negative, 
equal to for C > 1, where TZd is real and positive, and 
that is in between these two values for ( G [—1, 1], where 
TZd is complex. It follows from the symmetry of TZd that 



F{^0 = -n-F{C). 



(6.11) 



Figures || and ^ show F{Q for the square and cubic lat- 
tices. 



Us ing f (C), can be found from equations (4.17) 
and (F20I ) 



. int _ 
O2 — 



-4d/3t(C+l) 



(6.12) 



Partially integrating, we find 



tint 



Adf3t 



-Mpt 



8d(3t 



dCFiC) sinh 4df3tC 



(6.13) 



To find 62"*, the last equation can be numerically inte- 
grated for various values of f3t. 

As before , we compare this with the result for the ideal 
lattice gas ( 3.15| ) 



1 



bl = ^61(2/?) = -e-4'i/3t[/o(4/3t)]' 



(6.14) 



The ratio 62"*/fc§ is given by 

-2 



1 - ^ f dCFiO sinh 4d(3tC } . 



(6.15) 



For pt ^ 0, this ratio goes to —2. The second virial coef- 
ficient 62 = 62+^2"* thus approaches —62 for high temper- 
atures. This is indicative of fermionic behaviour, which 
is indeed what one would expect: at high temperatures, 
the only important contribution to the free energy of the 
system is the entropy involved in distributing a certain 
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number of hard core particles over the lattice. This is the 
same as for fermions. For fit ^ oo the behaviour of ( |6.15| ) 
depends on the behaviour of F{C,) for C ^ 1. In two di- 
mensions, oc (1 — C)/ ln(l — C), which leads to b™'^ 
going to zero like —1/ ln(/3i) as (it ^ oo. The same low- 
temperature behaviour is found for quantum hard disks 
In three dimensions, F(C) oc (1 — C)^, which gives 
1^2^ oc —1/y/pt for jit oo, just as for quantum hard 
spheres j^]. So in both cases, 62 approaches the value for 
the ideal Bose gas at low temperatures. This shows that 
when the thermal wavelength ( |3.9[ ) exceeds the lattice 
constant a{— 1), the effects of the Bose statistics starts 
to dominate. Plots of 62/62 versus fit for two and three 
dimensions are given in figures I and |. 



VII. CONCLUSION 

The pressure and density of a gas of quantum particles 
on a lattice can be expressed in fugacity expansions, 

^ " i=i 



n 



^lyebee-'^^. (7.2) 



For bosons with a hard-core interaction, the first two 
coefficients, 61 and 62, can be calculated. It turns out 
that the effect of the hard-core interaction on 62 depends 
strongly on temperature and on the transfer integral t. 
For small pt it is fermionic in character, as far as 62 
is concerned. For large (it the effect of the hard-core 
disappears, and only the Bose character remains. 

The system of vacancies that exists in solid ^He should 
behave to a good approximation like this simple model: 
the vacancies move through the crystal lattice by a tun- 
neling process, they are bosons, and they have a hard- 
core repulsion. The gas of vacancies has been experimen- 
tally studied by probing the attenuation and the velocity 
shift of sound in a ^He crystal Q. Very pure hep He- 
lium was used, so that the effects of both the phonons 
and the delocalized, bosonic vacancies could be observed. 
Treating the vacancies as a gas of free particles, it was 
seen that they obey Bose statistics. For the expressions 



(7.1) and (7.2) this means that not only the first term. 



but at least also the second is experiment ally o bservable 



One can use the expressions ( p.q ) and ( 3.15 ) for a gas 



of free, ideal bosons to e stim ate th e or der of magnitude 



of the various terms in ( [7.l| ) and (7^). One finds that 
6? = vo/l'^/^X^, so that 



^ 00 



/3A3 ^ £5/2 

^ e=i 



(7.3) 



1 e-'^f^^ 



X3A^ £3/2 
1=1 
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Using the value A/ks = 0.71K found in and the max- 
imum temperature T = 0.85K at which the experiments 
were performed, this sho ws t hat the ratio between the 



second and first terms in (7.4) is 0.15, and that between 



the third and first terms is 0.04. 

If experiments can be done tha t de tect the contribu- 
tion of 15% of the second term in ( [7.4[ ), it is also possible 
to detect the effects of the hard-core on the coefficient 62 



in (7.2). As can be seen from figure |6|, it varies consider- 
ably with pt, and it can even change sign. However, it is 
difficult to extract this information from the attenuation 
experiment, since there it is not clear what exactly the re- 
lation between the measured quantities and the vacancy 
density is. It would be necessary to directly measure, say, 
the pressure due to the vacancies S, and then compare 



this with (7.1). In such an experiment it would be cru- 
cial to take the hard-core effects into account. Work is in 
progress to calculate the coefficient 62 for the hep lattice; 
however, the results are not expected to differ much from 
those given for the simpler cubic lattice here. 

The results for the two-dimensional system might have 
relevance for low-density ^He films on textured substrates 
1^. There, the presence of adsorption sites localizes the 
Helium particles on the points of a two-dimensional lat- 
tice, defined by the substrate. Thus the tight-binding 
approach used in this paper becomes applicable. The 
adsorption sites cannot hold more than one particle, so 
that the hard-core repulsion is also present. 
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The integral (6.5) can be expressed in terms of elliptic 
integrals as follows ||l^,|ll[. TZ2 is equal to 

7e.(C,i?)~/ dk,^^ 

1 



I dky ^ , (5) 

where ky S [0, tt]. On making the substitution p = cos ky, 
and writing 

A^-B^ = '-^^^^i,.~p)i,.^p), (6) 
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where 



^ -2C + cos(jr J2) ^ -2C - cos(g,/2) 

cos(ifj,/2) ' cos(ifj,/2) ' ^ ^ 

we end up with integrals of the type 

dp^, (8) 



where w = {1 ~p){l +p){qi ^ p){q2 ~ p) and p G [—1, 1]. 
By defining 

q- = mm{qi,q2), q+ = ma.x{qi , q2) , (9) 

we see that values of p € g+J give A"^ < J3^, and thus 
contribute to the imaginary part of 7?.2. Values of p out- 
side this interval give > and thus contribute to the 
real part, for p < q- with sgn(A) = — sgn(cos(iiry/2)), 
and for p > q+ with sgn(yl) = sgn(cos{_i4ry/2)). The re- 
sult is a sum of integrals of the type (§) between limits 
that are consecutive zeroes of w. These integrals can all 
be expressed in terms of the complete elliptic integral of 
the first kind, K{x) pi, 



K{x)= dp , (.Te[0,i]). 



(10) 



There are several cases to be considered: 



1- '7-,g+ < -1 (la) or > 1 (lb) 

^ At^CO^{Ky/2) V(q+-l)(q_ + l)' ^ ^ 
where e = + for (la) and e — — for (lb). 



2. g_ < -l,g+ > 1 or G [-1,1] 

^'(^'^^ = cos{Ky/2)\ y(^l)(l - ^^^^ 



3. q- < -l,g+ G [-1,1] (3a) 
or g_ G [-1,1], g+ > 1 (3b) 



4i^cos(i^y/2) ^2{q+ - q_) 
I K{l/r2) 



4t7rlcos(if,/2)l ^2(g+ -q_) 



(13) 



where e = -f for (3a) and e = — for (3b). In the above, 
ri and r2 are given by 
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r2 



{q+ + m-< 



(14) 
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FIG. 1. The Brillouin zone for the one-dimensional prob- 
lem. Domains I and HI' are the Brillouin zone used in the 
calculation. 



FIG. 2. The phase #(C, K) in the C, K/2 plane. 



FIG. 3. The function F{() for the square lattice. The only 
values of interest are those between ^ = — 1 and C = Ij where 
F{C,) is not constant and varies from — tt to 0. Only positive 
values of ( are shown, since -F(C) for negative ( can be found 
using the antisymmetry of that function around C = and 
F(C) = -7r/2. 



FIG. 4. The function F{Q for the cubic lattice. 
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FIG. 5. The ratio 62/62 as a function of /3t for the square 
lattice. For high temperatures (small /3t), 62 = —62; which 
is the same value as for formions. For lower temperatures its 
value crosses over to 62 = 62; the ideal boson value. 

FIG. 6. The ratio 62/62 a-s a function of pt for the cubic 
lattice. Note that the cross-over to the ideal boson value is 
much faster than for the square lattice. 
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